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PACS 03.65.Vf - Phases: geometric; dynamic or topological 
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Abstract. - We study the transport of cold fermionic atoms trapped in optical lattices in the 

presence of artificial Abelian or non-Abelian gauge potentials. Such external potentials can be 
created in optical lattices in which atom tunneling is laser assisted and described by commutative 
or non-commutative tunneling operators. We show that the Hall-like transverse conductivity of 
such systems is quantized by relating the transverse conductivity to topological invariants known 
as Chern numbers. We show that this quantization is robust in non-Abelian potentials. The 
different integer values of this conductivity are explicitly computed for a specific non-Abelian 
system which leads to a fractal phase diagram. 



During the last decades, remarkable riuantum phenom- 
ena have been discovered such as the Aharonov-Bohm ef- 
fect [1], the geometric phases [2,3], and the integer quan- 
tum Hall effect (IQHE) [4, .5]. These important phenom- 
ena manifest themselves when gauge potentials are present 
in non-trivial topological spaces, in which cases the wave 
functions describing the system may acquire topological 
or geometrical phases. They are elegantly described in 
terms of differential geometry and topology in which the 
topological phases are related to the concept of holon- 
omy [6]. In the IQHE, Hall's transverse conductivity 
of two-dimensional electronic systems submitted to high 
magnetic fields turns out to be quantized according to 
(Jh = —C (e^//i) where C is a topologically invariant in- 
teger called the first Chern number [6-8]. Since the pro- 
portionality factor is a combination of the electric charge 
e and Planck's constant /i, such quantization phenomena 
play a key role in metrology [4]. 

Very recently, great experimental advances have been 

performed on the control of cold atomic systems [9-12]. 
In such systems, artificial gauge potentials can be created 
in optical lattices in order to reproduce the dynamics of 
periodically constrained fermionic atoms submitted to the 
analogue of a magnetic field. Such systems would thus of- 
fer the possibility to observe the corresponding fractal en- 
ergy spectrum known as the Hofstadter butterfly [13, 14]. 
More recently the realization of non-Abelian gauge po- 



tentials has also been envisaged, allowing the observa- 
tion of a non-Abelian Aharonov-Bohm effect [15], mag- 
netic monopoles [16], and particular metal-insulator tran- 
sitions [17]. 

The purpose of this Letter is to report the possibility 
of an integer quantum Hall-like effect in the transport of 
cold fermionic atoms trapped in optical lattices submit- 
ted to artificial gauge potentials. We consider this effect 
in a family of non-Abelian gauge potentials, which con- 
tains Abelian potentials as particular cases. In spite of 
the new topological context provided by the non-Abelian 
gauge structure, we show that an integer quantum Hall- 
like effect is possible and that the transverse conductivity 
is indeed quantized in terms of Chern numbers in sys- 
tems of fermionic cold atoms in optical lattices with ar- 
tificial Abelian or non-Abelian gauge potentials. Since 
the current induced in our system is electrically neutral, 
this quantized conductivity has the distinctive units of the 
inverse of Planck's constant, which could be important 
for metrological purposes. In this Letter, we first define 
the system Hamiltonian characterized by its non-Abelian 
gauge structure and evaluate the transverse conductivity 
with Kubo's formula which gives the linear response of 
the system to an external static perturbation. We then 
express the Hall-like conductivity in terms of Chern num- 
bers in order to prove the quantization of this physical 
quantity. We compute these topological invariants for a 
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specific system which leads to a fractal phase diagram. We 
then disc;uss possible experimental setups for the valida- 
tion of our theoretical considerations. 

We are interested in systems described by the following 
Hamiltonian 

} + h.c. 

(1) 

which determines the dynamics of fermions in a two- 
dimensional lattice within a tight-binding approximation 
as studied experimentally in cold atom systems trapped in 
optical lattices [9-12]. For Ux,y belonging to the Abelian 
group of unitary complex numbers U{1), this Hamilto- 
nian reproduces the evolution of an electronic system con- 
strained by a periodic potential and submitted to a mag- 
netic field [14]. In the latter case, the single-particle wave 
function gets a phase factor Ux {Uy) when tunneling is per- 
formed along the x (y) (iirec;tion and their pro(luc;t around 
a unit cell gives a total phase proportional to the pene- 
trating magnetic flux. In this Letter, we study the case in 
which the system has a U(2) gauge structure, which im- 
plies that the operators Ux and Uy are 2x2 unitary matri- 
ces and the single-particle wave function has two compo- 
nents: when tunneling occurs along the x (y) direction, the 
matrix Ux {Uy) acts as a tunneling operator on the two- 
component wave function and the system is non-Abelian 
for [Ux,Uy] 7^ [15,17]. The operators U^.y are related to 
the gauge potential A present in the system via the usual 
relation Ux^y = e*"^^'". In order to indicate multicompo- 
nent wave functions and operators, we use the following 
notations introduced in Ref. [3]: A row of n orthogonal 
kets is denoted by {[V")! = •••7 IV')n) and a column of 

n kets by ||'0)}.The coefficients ta and tb describe trans- 
fers in the x = mix and y = nly directions, where {Ix, ly) 
are the unit cell's lengths, and |cj„ „} (resp. |cm,ri}) are 
the 2-component fermionic creation (resp. annihilation) 
field operators on site (m,n). After a change of basis, 
the latter's components are transformed into Cm,nj = 
X^A VxjijnTn) c\, where c\ is the operator of annihilation 
of a fermion in the j^^ component tpxj{in, n) of the single- 
particle wave function {'tp\{r)\ = (V'Ai(^); V'A2('")) with 
r = (m, n) and j = 1,2. 

The current density of the system „ associated with 
the Hamiltonian operator (1) is written in terms of the 
fermionic operators and its components are given by 



h 



J2 ^A,(™ 

AA'jj' 
+ h.C. 



^,n){Ux) ..,tp\'j'{m,n)c\cx' 



and a similar expression for jm,n,y 

It has been recently suggested that non-Abelian gauge 
potentials can be created in optical lattices [15, 17]. In 
such setups, atoms with doubly degenerate Zeeman sub- 
levels hop from a site to another with the assistance of ad- 
ditional lasers and their tunnelings are indeed described 



by unitary non-commutative operators [15]. The corre- 
sponding gauge potential is given by 



A = 



' 1 2^ 

— i(f> TT 



/ 27rmai 

'V 





27rma2 



(2) 



and induces an effective "magnetic" field characterized by 
three parameters ai, a2 and (j), which can be controlled 
with external lasers [15]. We notice that this gauge po- 
tential is Abelian under the condition that ai — 0.2 is an 
integer and non-Abelian otherwise. 

We evaluate the linear response of the system to an ex- 
ternal static force applied along the y direction Tiext = 



}. A generalized Lorentz force 
produced by the effective "magnetic" field induces a trans- 
verse current which satisfies (ja,) = Oxy Ay. Kubo's for- 
mula expresses the transverse conductivity in terms of the 
total current J = Em,n Jm,n 



^ ih^ {h\Jx\'S){^\Jy\h) 



{0\Jx\b){b\Jy\0) 

{Eb - EqY 



(3) 



where \h) denotes the eigenstates of the Hamiltonian (1) 

of eigenvalue Eb and V is the system volume or hyvh. 

The system presents symmetries under translations de- 
fined by the operators \{T^^){m,n)} = |?/'(m -I- g,n)} 
and \{Ty^){m,n)} = |'0(m, n+ 1)}. The operator Ty 
commutes with the single-particle Hamiltonian H, since 
the potential (2) only depends on the x coordinate. On 
the other hand, the commutation relations \T^,H\ = 
\T^,Ty\ = are satisfied under the conditions aj = pj/q 
with j ~ 1,2 and pi,p2,q some integer numbers. The 
translational symmetries allow us to write [■i/'(TO,n)} = 
e*'^''|w(m)} with \u{m)} g-periodic and the first Bril- 
louin zone is a 2-torus defined by fc^, e [0, ^] and 
ky G [— 7r,7r]. The latter change can be viewed as a gauge 
transformation which modifies the Hamiltonian by substi- 
tuting An ^ A'^ = An + kn so that the single-particle 
Hamiltonian takes the following form 

H{k)\u{m)} = 

i^ei(A.+fe»)|^(^ + 1)} + i„e-'(^-+*=»V(m - 1)} 
+ tbe'^^y+''y^\uim)} + tbe-'^^y+^y^ \u{m)) (4) 

This equation can be differentiated with respect to k^ and 
a straightforward calculation yields 

^ {Eb-E,y -¥ ^ ^ {{uy\dk^H\u,)} 



eA<eF e;^/>€F 



x{{ux\dk^H\uy)}/{ex-eyf (5) 

where H\ux} = £x\ux} are the single-particle eigen- 
states and the scalar product is defined by {(wa|wa')} = 

I]m=iI]j=i,2"Aj(™)"A'j('7i). The Fermi energy ep is 
supposed to be situated inside a gap of the spectrum. We 
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find that axy — —(^yx and the identity (5) allows us to 
write the conductivity as 



'xy 



(2 



dk 



(6) 



This result is important because it relates the transverse 

conductivity to integer topological invariants called Chern 
numbers. The latter are defined on a fibre bundle which is 
a topological space that locally resemble the direct prod- 
uct of the parameter space with the non-Abelian gauge 
group U{2) and that we note P{T^,U{2)). This prod- 
uct space may be globally non-trivial and the structure 
twists when the Chern number is non-zero. In the Abelian 
framework, the interpretation of the IQHE in terms of 
topological arguments [6,7] showed that Berry's curvature 
= [{dkMk)\dkM'')) - {dkjmdkJKk))]dk,dky, 
with \tp) an eigenstate depending on the wave vector fe, 
defines a curvature on the fibre bundle P{T^,U{1)) and 
that its integral over the parameter space corresponds 
to the first Chern number C = ^ Jj2 ^ ^ which is an in- 
teger. These concepts can be generalized in the actual 
non-Abelian framework [18, 19] for which Berry's curva- 
ture is a 2 X 2 matrix and the Chern number is then 
defined by 

c=^ J ^tr{\{dk^u\}{\dk,u)\ - \{dk,u\}{\dk^u)\)dk,dky 

(7) 

where ti^F denotes the trace of the matrix J^. One even- 
tually finds that the Hall-like conductivity (6) is given by 
a sum of Chern numbers 



(8) 



eA<eF 



where the integer numbers C{ex) are associated with each 
energy band. It follows that the transverse Hall-like con- 
ductivity of the system evolves by steps corresponding to 
integer multiples of the inverse of Planck's constant and 
that it remains constant under small perturbations. Be- 
sides, a measure of the transverse conductivity in this non- 
Abelian system provides a direct evaluation of these topo- 
logical invariants. 

The Chern numbers were explicitly obtained in the 
usual Abelian framework [20], and we here show that they 
can be computed explicitly in non-Abelian systems as well. 
Using Stokes' theorem one gets the following result 



^ xy 



(9) 



where 6\ is the phase accumulated by u\ around the torus 

and evaluated on its border 9T^. These phases can be 
evaluated in a way similar as in the Abelian case [20]. 



Let us write the Schrodinger equation associated with the 
non-Abelian Hamiltonian (4) in the form 



ta 



+ tb 



^'•'U2{m — 1) 
-i<t>+ik.ui{m + e-''f'-^^-ui{m - 1) 

2cos{2TTaxm + ky)u\{m)\ _^(ui{m) 



2 cos(27rQ;2'TT' + ky)u2{m) 



U2{m) 



(10) 



When fa — » 0, this equation gives a 2g-degenerate single 
band and the dispersion law for Uj{m) is 

ej{m,ky) = 2tb cos{2Trp jm/q + ky) =: 2t(,cos(A;™) (11) 

First, we consider the case where ai = pi/q and 
C(2 — {pi + nq)/q with n integer. This case allows the 
explicit computation of Chern numbers for Fermi energies 
situated inside the various gaps of the spectrum in the ap- 
proximation of weak coupling with ta small and tb = 1- 
Since the Chern numbers are topological invariants, their 
weak-coupling values extend to other couplings as long as 
the band gaps deform continuously. As ta increases from 
to 1, different regimes are observed. For > and 
small, q — 1 gaps open at the crossing points ky = ±y 
with I integer. Thereafter, the dispersion branches sepa- 
rate and no crossing remains. When ta = I, the spectra 
e = e{ai = ^;a2 = 2ii^^) precisely depict a Hofstadter 
butterfly [13]. Finally when ta > I, q new gaps open so 
that the number of gaps in the spectrum becomes 2q — 1. 

In order to obtain the Chern numbers, one can study 
how the eigenfunctions change at the different gaps as ky 
varies in the flrst Brillouin zone [20]. Near these gaps, 
the four functions u,i.2(w.) and '(ii,2("''') with their disper- 
sion branch crossing each other in the limit ta = are 
strongly coupled together. Anticrossing appear for ta ^ 
which can be calculated by perturbation theory giving an 
effective Schrodinger equation for the vector v{m,m' ,r) 
formed with these four functions. Considering the r**^ gap 
and setting tr = m — m' , this equation can be written as 
M{q, tr, kx) v{m,m',r) = €v{m, m',r). The latter system 
has been solved in order to find the four eigenvectors Vi. 
When ky passes a crossing point, the eigenvectors undergo 
some transformation Vi Vi = R{i,tr, kx)vi. Following 
the eigenvector of the r**^ band through the successive an- 
ticrossings in the first Brillouin zone ky e [—n,n], one 
finds that the eigenvector goes under the total transfor- 
mation \u\} — » e*^*-*''"*''"^-*'^"' |7iA}- We can now evaluate 
the r**^ band's contribution to the transverse transport 
coefficient axy by considering the formula (9) which gives 
(Txyir^^hanad) = {tr — tr-i)/h. 

The integers tr have still to be evaluated. From the def- 
inition (11), one has ei{m,,ky) = e2{m,ky) in the present 
situation where a2 = {pi + nq)/q. Two conditions are 
necessary for crossing points that determine the r*^'^ gap 
when ta 0: fc™ = —k"^ + 2TTSr with Sr integer, and 
fc^' = fc^ — 2'Kpitr/q with tr = m — m', which lead to the 
Diophantine equation r = pitr + qSr, with r the position 
of the gap and ai = pxj q,oi2 = {pi+nq)/q- For ta > 1, all 
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Fig. 1: Energy spectrum e — e(ai) for 02 = 31/701. Gaps are 
shown in black. A Hofstadter butterfly crosses this structure 
at Qi = a2 = 31/701, where its larger gaps are put forward 
through thin black lines. 

the 2q — l gaps open, so that axy = X]5=i(*i ~ ^j-i)/'^ = 
tr/h, if the Fermi energy lies in the r^^ gap. For ta < I, 
only q ~ I gaps open and the bands have a doubled mul- 
tiplicity. Accordingly, one finds axy — 2 tr/h if the Fermi 
energy lies in the r"^ gap. 

When resolving the non-Abelian Harper equation (10) 
{ta = tb = 1) one finds the energy spectrum e = e(ai,a2) 
which forms, in the 31? space (e, ai, (22), a structure called 
the Hofstadter moth [15]. This 3D structure is character- 
ized by a complex distribution of small gaps as shown in 
Fig. 1 which illustrates a typical cut through the "moth" . 
If the spectrum e = e(ai) is depicted in a particular plane 
a2 = (pi + nq) I q with n integer, one recovers the Hofs- 
tadter butterfly [13,21-23]. In this spectrum, each gap is 
characterized by the two integers Sr and tr given by the 
aforementioned Diophantine equation [23]. Consequently, 
it is possible to draw a phase diagram similar to Osadchy 
and Avron's [24] for this new non-Abelian case. This new 
phase diagram represents the integer values of the trans- 
verse transport coefficient (Jxy inside the infinitely many 
gaps of the butterfly, which are given by Oxy ~ 2 t^jh. In 
order to check the robustness of this result, we depict in 
Fig. 2 the energy spectra for different values of a\ — a2 and 
anisotropy ratio t^jta- The gauge potential is Abelian in 
Figs. 2(a) and (b), but non-Abelian in Figs. 2(c) and (d). 
We observe in this latter case that the spectrum thickens 
in particular around a\ — 0.5. Under the effect of the 
perturbations, some gaps disappear as in the Hofstadter 
model [25] , but the principal gaps where the Hall- like con- 
ductivity is quantized remain. 

An important issue is whether this result could be ex- 
perimentally probed. A possible system described by the 
Hamiltonian (1) which includes the gauge potential (2) 
has already been proposed by Osterloh et al. [15] who put 
forward the possibility to submit cold fcrmionic atoms to 



such artificial non-Abelian gauge potentials. Each vertex 
of a 3D optical lattice contains an atom with degenerate 
Zceman sublevels in the hyperfinc ground state manifolds 
{\g)j, \e)j} with j — 1,2, so that the states we are deal- 
ing with are represented by the rows {|e)| = (|e)i,|e)2) 
oi' {Iff) I = ( Iff) 1 J Iff) 2)- We note that fermionic atoms 
with such properties, i.e. '^^K in states \F = 9/2, mp = 
9/2, 7/2, ...) and \F = 7/2, mp = -7/2, -5/2, ...) as sug- 
gested in Ref. [15], are already studied experimentally in 
optical lattices [10-12]. Lasers are used in order to cre- 
ate the lattice, as well as the non-Abelian gauge potential 
through state-dependent control of hoppings which take 
place within every plane z = constant. In this frame- 
work, we consider a possible setup that exhibits a cur- 
rent in the y direction consisting of an accelerating optical 
lattice [14,26], with the (possibly gravitational) accelera- 
tion a ~ {0,ay). Consequently, a perturbative term Hcxt 
with the intersite energy difference Aj, = Mayly where 
M is the mass of the particles is added to the Hamilto- 
nian [14,27]. In this case, the transverse transport coeffi- 
cient fjxy would give the relation between the external forc- 
ing and the transverse atomic current through the lattice. 
We notice that the current is electrically neutral so that 
this conductivity has the distinctive units of the inverse of 
Planck's constant and relates the energy difference Ay to 
a number of atoms per unit time in the current. Recently, 
it has been suggested that such atomic currents could be 
experimentally produced and measured in optical lattices. 
In particular, Ponomarev et al. [27] have investigated the- 
oretically the possibility of such currents in a system of 
fermions coupled to bosons. The latter acts as a cold bath 
for the fermions, which induces a relaxation of the Bloch 
oscillations. If a transverse current is present in such sys- 
tems under the conditions described in this work, we have 
shown that the transverse Hall-like conductivity should 
be quantized. In order to observe experimentally the ef- 
fects emphasized in this work, it is crucial that the energy 
gaps represented in Fig. 2 remain sufficiently clear under 
thermal effects. The energy resolution required is approx- 
imatively 100 Hz which corresponds to temperatures of 
the order of 10 nK. These temperatures are realized in ex- 
periments involving '^^ K in state \F = 9/2, mp = 9/2) as 
mentionned in Refs. [10-12]. We finally point out that re- 
cent studies put forward the possibility to create and mea- 
sure currents in atomic periodic systems, such as optical 
lattices, by connecting two reservoirs in order to create a 
chemical potential gradient through the system which gen- 
erates a current [28]. Such devices, in which the trapping 
potential present in the optical lattice can be taylored, are 
envisaged in the new context of atomtronics [28] and could 
be considered for the detection of the quantized Hall-like 
conductivity in optical lattices submitted to Abelian or 
non-Abelian gauge potentials. 

In this Letter, we have shown the existence of topologi- 
cal invariants in the transverse transport properties of cold 
atomic systems. It appears that fermionic atoms submit- 
ted to artificial Abelian or non-Abelian gauge potentials 
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Fig. 2: Spectra for ai — pi/701 in the Abelian cases: (a) 02 = pi/701 + n with A — 2tb/ta ~ 2; (b) Q2 — pi/701 + n with 
A = 3.5; and non-Abehan cases: (c) 0:2 = {pi + 0.1)/701 + n with A = 2; (d) 02 = {pi + 0.05)/701 + n with A = 2.5. The 
integers in the energy gaps give the corresponding values of hay^- 
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and perturbed by an external forcing should exhibit an in- 
teger quantum Hall-like effec;t. We have here shown that 
the quantization of the transverse conductivity is robust 
if the gauge potential is set non-Abelian. The different 
values of the Hall-like conductivity are computed for a 
specific U (2) system which depict a fractal phase diagram 
that should partly resist perturbations encountered in a 
realistic setup. Cold atom systems trapped in optical lat- 
tices present the versatility propitious to the study of the 
fractional quantum Hall effect (FQHE) [29]. These sys- 
tems also constitute good candidates for the exploration 
of the quantization and fractal structures we have here 
emphasized. We can speculate that the quantization of 
the transverse conductivity of neutral fermionic cold atom 
currents could open a new window in metrology. 
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